Comparison of Energy Efficiency between Atmospheric Batch Pressure-Retarded Osmosis and Single-Stage Pressure-Retarded Osmosis

Batch pressure-retarded osmosis (PRO) with varied-pressure and multiple-cycle operation using a pressurized variable-volume tank has been proposed as a high-efficiency osmotic energy harvesting technology, but it suffers scalability constraints. In this study, a more scalable batch PRO, namely, atmospheric batch PRO (AB-PRO), was proposed, utilizing an atmospheric tank to receive and store the intermediate diluted draw solution (DS) and a pressure exchanger to recover the pressure energy from the diluted DS before being recycled into the tank. Its performance was further compared with single-stage PRO (SS-PRO) at different flow schemes via analytic models. The results show that the AB-PRO with an infinitesimal per-cycle water recovery (r) approaches the thermodynamic maximum energy production under ideal conditions, outperforming the SS-PRO with lower efficiencies caused by under-pressurization (UP). However, when considering inefficiencies, a ~40% efficiency reduction was observed in AB-PRO owing to UP and entropy generation as the optimal r is no-longer infinitesimal. Nonetheless, AB-PRO is still significantly superior to SS-PRO at low water recoveries (R) and maintains a stable energy efficiency at various R, which is conducive to meeting the fluctuating demand in practice by flexibly adjusting R. Further mitigating pressure losses and deficiencies of energy recovery devices can significantly improve AB-PRO performance.


Introduction
Renewable energy, as an alternative to fossil fuels in power production to reduce carbon emissions, has gained considerable attention in recent decades [1,2]. Among various types of renewable energy, osmotic energy, also called salinity-gradient energy, which originates from the mixing of two solutions with different salinities, has raised interest owing to its huge global capacity and potential accessibility from both natural and industrial sources [3][4][5][6]. The global osmotic power generated from the mixing of river water with seawater was estimated to be 1.6-2.6 TW, of which around 980 GW could be extracted as electricity, which is comparable to that generated from hydropower [7]. Despite the enormous potential, osmotic energy has not been tapped at large scales.
Among various technologies for the harvesting of osmotic energy, membrane-based pressure-retarded osmosis (PRO) [8,9] and reverse electrodialysis (RED) [10,11] have the highest potential for large-scale applications. Simulated and experimental results from previous studies indicate that PRO outperforms RED in terms of higher membrane power density [12][13][14]. In PRO, a hydraulic pressure, lower than the osmotic pressure difference In this study, an alternative design of the BPRO by incorporating atmospheric tanks for storing the intermediate DS and FS and a pressure exchanger for recovering pressure energy from the DS effluent as shown in Figure 1c, namely, atmospheric batch PRO (AB-PRO), was proposed, which avoids the use of intricate pressurized variable-volume draw tank in the BPRO reported previously. Specific energy production (SEP) and energy production efficiency (EPE) of AB-PRO were calculated and compared with other PRO device (i.e., PX or hydro-turbine), π D,in and π D,out : osmotic pressure of DS influents and effluents, respectively, π F,in and π F,out : osmotic pressure of FS influents and effluents, respectively.

Configuration and Operation
As shown in Figure 1a,b, SS-PRO is operated in an open-loop, continuous mode, where both the DS and the FS are discharged out of the system after passing through the membrane module once. The DS influent is pressurized through exchanging hydraulic pressure with the DS effluent inside a pressure exchanger (PX). Owing to the imperfect efficiency of the PX, a circulation (or booster) pump is installed to further pressurize the DS influent to the target value before it enters the membrane module. Under a constant hydraulic pressure below the osmotic pressure difference (∆π), water molecules permeate through the membrane from the FS to the DS, which results in a diluted DS effluent with an increased flow rate. The diluted DS effluent is split into two streams. The first stream, at the same flow rate as the DS influent, is directed to the PX to exchange pressure with the DS influent flowing to the module. The other stream at the same flow rate as the water permeating through the membrane is depressurized by an energy recover device (ERD), such as a hydro-turbine (HT) or PX, to produce electricity or pressurize a fluid in anther system [28][29][30]. The DS and the FS flowing to the membrane module can follow two different flow schemes, the co-current mode ( Figure 1a) and the counter-current mode ( Figure 1b).
In contrast, as shown in Figure 1c, AB-PRO is operated in a closed-loop, variablepressure mode with two atmospheric tanks to store DS and FS. In this case, the intermediate diluted DS and the concentrated FS are recycled and mixed with the residual solutions in the DS and FS tanks, respectively, instead of being discharged out of the system. It allows AB-PRO to be operated in multiple cycles with a smaller water recovery in each cycle (i.e., per-cycle water recovery) to achieve the total water recovery requirement. Therefore, the AB-PRO process can start at a higher hydraulic pressure because only a small portion of water passing through the membrane from FS to DS in the first cycle brings about a milder decline in the osmotic pressure difference along the membrane module than that of the SS-PRO, which is also demonstrated by the osmotic pressure profiles (dotted curves) in Figure 1. As DS is diluted and FS is concentrated from cycle to cycle, the applied hydraulic pressure in AB-PRO is gradually reduced to the same level as the constant pressure applied in SS-PRO. Owing to a higher average applied pressure, AB-PRO features less energy loss due to under-pressurization compared to SS-PRO. Moreover, the AB-PRO proposed in this study adopts a PX for the energy exchange between the influents and effluents of DS, which allows the practical atmospheric tank to be used instead of the less scalable pressurized variable-volume tank proposed in the previous study. Both co-current and counter-current flow orientations can be employed in AB-PRO, but the impact of flow schemes is insignificant since the change of osmotic pressure difference in each cycle of AB-PRO is not obvious due to the small per-cycle water recovery (r) as illustrated by the osmotic pressure profile in Figure 1c.

Derivation of Energy Production
To have an insight into the practical performance of each PRO process, the impacts of inefficiencies of devices (e.g., pump, PX, and ERD) and pressure losses on both the DS and FS sides were systematically assessed. Analytic expressions of the practical maximum specific energy production (SEP) and energy production efficiency (EPE) were derived to evaluate the energy production performance of different PRO modes under ideal and practical conditions. SEP is defined as the energy generated per total volume of DS and FS, and EPE refers to the ratio of the energy extracted by PRO to the Gibbs free energy generated from the mixing of DS and FS. To simplify these analyses, it was assumed that the salt rejection of the membranes is 100%. The osmotic pressure (π) was assumed to be linearly proportional to the salt concentration (c), which can be expressed by the van't Hoff equation [27,31]: where v is the number of ionic species each solute molecule dissociates into, R g is the ideal gas constant, and T is the absolute temperature, which is assumed to be 298 K in the following calculations.

Specific Gibbs Free Energy
The Gibbs free energy of mixing (∆G m ) normalized by the total volume of DS and FS defines the thermodynamic maximum specific energy production (SEP theomo ). As illustrated in Equations (2)-(4) [19,32], the SEP theore is related to the concentration of FS and DS, the water recovery, and the volumetric fraction of DS and FS.
where V D,0 and V F,0 are the initial volumes of DS and FS, respectively, π D,0 and π F,0 are the initial osmotic pressures of DS and FS, respectively, R is the water recovery defined as the ratio of total water permeation volume (V P,T ) to V F,0 , which ranges from 0 to 1 (when FS is pure water), and ∅ is the volumetric fraction of FS. It should be noted that ∅ was optimized according to Equation (A1) in Appendix A in the following calculations.

Energy Production Performance of SS-PRO
In SS-PRO, the highest constant hydraulic pressure (P C ) equals to the lowest osmotic pressure difference (∆π) along the membrane module. For a given R and ∅ in SS-PRO, the practical maximum SEP (SEP SS−PRO ) is obtained at the highest applied constant pressure. Therefore, SEP SS−PRO can be computed via Equation (5), and the SEP under ideal conditions (SEP SS−PRO,ideal ) can be obtained by Equation (6).
where η P , η PX , and η ERD are the efficiencies of pump PX and ERD, respectively, and P D,loss and P F,loss are the pressure losses on the DS side and the FS side, respectively. Equation (5) also reveals the four energy components contributing to SEP, the specific energy production through the ERD (η ERD RP C ), the specific energy consumption for pumping (− η P P C ), the specific energy recovered by the PX ( , and the specific energy loss due to pressure drop ( ∅ η P P F,loss + (1−∅) η P P D,loss ). P C is different with various flow schemes in SS-PRO. As shown in Figure 1a, the highest pressure in the co-current SS-PRO (P C,co ) depends on the ∆π between the final diluted DS and the final concentrated FS at the end point of the membrane module (Equation (7)): However, the highest pressure in the counter-current SS-PRO (P C,counter ) is determined by the smaller one between the ∆π at the two ends of the membrane module (i.e., at the outlet of FS or DS).
(1−R) , at the outlet o f FS and when P 1 < P 2 − π F,0 , at the outlet o f DS and when P 1 > P 2 (8)

Energy Production Performance of AB-PRO
Unlike the constant pressure applied in SS-PRO, AB-PRO features variable-pressure according to the changing ∆π during the recirculation of DS and FS. Therefore, the hydraulic pressure applied in AB-PRO is primarily related to the per-cycle water recovery (r), rather than the total target water recovery (R). In this study, r was assumed to be constant during a batch of AB-PRO. Moreover, the spatial effects in AB-PRO were ignored by assuming a spatially invariant concentration throughout the pressure vessels and tanks at a given moment. Such an assumption facilitates the derivation of the analytic expression of SEP of the AB-PRO without compromised accuracy since the length of the pressure vessel in AB-PRO can be very short and the change of ∆π can be finished in a short time [27]. The impact of flow schemes on SEP was also ignored based on this assumption. The maximum variable pressure applied in AB-PRO is a function of the time (t) as follows: where Q F and Q D are the FS and DS circulation rates, respectively, Q P is the water permeation rate, τ is defined as the maximum retention time of the FS under the permeate flow rate of Q P , θ is the ratio of Q F to Q D in AB-PRO, which was assumed to be the same as that in SS-PRO. The approximation of practical maximum SEP of AB-PRO can be obtained via Equation (13). It should be noted that the pressure loss on both sides of DS and FS in AB-PRO is assumed to be 1/7 times that of SS-PRO due to the shorter module length of AB-PRO.
By differentiating the SEP AB−PRO with respect to r, the optimal r can be computed by the iterative method. Under ideal conditions, the maximum SEP (Equation (A14) in Appendix B) can be obtained with r → 0 , which equals the specific Gibbs free energy of mixing (Equation (2)).
The SEP of different PRO processes was compared under different scenarios. The values of parameters, including η P , η PX , η ERD , P D,loss , and P F,loss , are listed in Table 1. efficiency, P loss is the pressure loss. b the configured ERD can be either a hydro-turbine (HT) to generate electricity or a PX to exchange pressure with another system. c the pressure loss on both DS and FS sides is assumed to be the same (i.e., P D,loss = P F,loss ).  is the pump efficiency, is the pressure exchanger (PX) efficiency, is the energy recovery device (ERD) efficiency, is the pressure loss. b the configured ERD can be either a hydro-turbine (HT) to generate electricity or a PX to exchange pressure with another system. c the pressure loss on both DS and FS sides is assumed to be the same (i.e., , = , ). Figure 2 shows the SEP and EPE of different PRO processes with various pairs of DS and FS. The results indicate that the SEP increases with a higher salinity gradient between the DS and the FS in all three PRO processes.   figure (a-c), various salt solutions with different concentrations were used for osmotic energy harvesting including a 0.6 M NaCl solution (e.g., equivalent to the salinity of seawater) and a 1.2 M NaCl solution (e.g., equivalent to the salinity of SWRO brine) as draw solution (DS) alternatives, a 0.01 M NaCl solution (e.g., equivalent to the salinity of river water) and a 0.05 M NaCl solution (e.g., equivalent to the salinity of wastewater concentrate) as feed solution (FS) alternatives. It should be noted that the x axis, water recovery, in (a-c) is different from that in (d-f). (a-c) show the SEP and EPE performance at various target total water recoveries, while (d-f) indicate the change of osmotic pressure difference and hydraulic pressure at the real-time water recovery from 0 to the target value (R').

Ideal Energy Production Performance
As illustrated in Figure 2a,b, both co-current flow and counter-current flow SS-PRO processes show imperfect energy production efficiency (EPE) even under ideal conditions. The constant-pressure operating mode of SS-PRO results in unavoidable energy loss caused by under-pressurization (UP), which is indicated by the yellow areas in Figure 2d,e. The energy loss increases with the increase of R, leading to the decrease of EPE. However, there are still differences between the two SS-PRO processes, which have different flow schemes. For co-current SS-PRO, as R increases, the SEP first increases to a peak then decreases. According to Equation (6), without considering the inefficiencies, the SEP of SS-PRO only depends on (1) the applied pressure (P C ), (2) the water recovery (R), and (3) the volumetric fraction of FS (∅). For a specific pair of DS and FS at the optimal ∅ (Equation (A1) in Appendix A), when R increases, there is a tradeoff between the increased loss of energy owing to the drop of P C and the increased release of Gibbs free energy of mixing (∆G m ) due to more water permeation. When the adverse impact of UP dominates, both SEP and EPE drop significantly with the increasing R as a result of the declining P C owing to the rapid decrease of osmotic pressure difference (∆π) along the module in the co-current SS-PRO. When R is maximized, the osmotic pressure of the final concentrated FS is the same as that of the final diluted DS, which means that no hydraulic pressure can be applied to extract energy although ∆G m reaches maximum. Unlike the trend of SEP, the EPE keeps decreasing as R rises. Theoretically, to obtain a high SEP and EPE simultaneously, it is necessary to operate the co-current SS-PRO at low R using a pair of DS and FS with a high salinity gradient. For example, the SEP of 0.10 kWh·m −3 and the EPE of 65% at R = 0.50 using a 0.6 M NaCl solution as the DS and a 0.05 M NaCl solution as the FS can be improved to 0.14 kWh·m −3 and 86%, respectively, by changing the DS into a 1.2 M NaCl solution and reducing the R to 0.20. For the counter-current SS-PRO, the flow scheme induces a more even distribution of ∆π along the membrane module. As illustrated by the red dotted curve in Figure 2e describing the change of ∆π in the counter-current mode, it starts from a lower value and ends at a higher value compared to the black dash curve showing the change of ∆π in the co-current mode. It can be explained by the osmotic pressure difference profile depicted by the dotted curves in Figure 1c that the opposite directions of DS and FS flows in the counter-current SS-PRO result in a lower ∆π of (π D,in − π D,out ) at the DS inlet than that (∆π = π D,in − π F,in ) in the co-current SS-PRO, and a higher ∆π of (π D,out − π D,in ) at the DS outlet than that (∆π = π D,out − π F,out ) in the co-current SS-PRO. Since the hydraulic pressure applied in SS-PRO is constrained by the lowest ∆π, the counter-current flow mode contributes to a higher appliable hydraulic pressure and more extractable osmotic energy at the same R. Equations (7) and (8) also demonstrate a higher P C,counter than P C,co from a mathematical perspective. Therefore, the counter-current SS-PRO outperforms the co-current SS-PRO on both SEP and EPE especially at high R. In addition, as a considerable concentration change is required to bring ∆π to 0 (e.g., π D,in = π F,out , or π D,out = π F,in ) in the counter-current SS-PRO, the situation where no pressure can be applied to harvest energy will not occur in most cases. As shown in Figure 1b, the SEP of counter-current SS-PRO is monotonically increasing with R, and the EPE remains above 70% at any R.
In contrast, as shown in Figure 2c, AB-PRO can approach the thermodynamic maximum SEP and 100% of EPE under ideal conditions, which are significantly higher than those of the two SS-PRO processes. The mechanism illustrated in Figure 2f reveals that the energy loss caused by UP can be eliminated in the AB-PRO by operating with an infinitesimal per-cycle water recovery (r), where the variable hydraulic pressure applied on the DS side is as close as possible to the changing osmotic pressure difference. Using the 1.2 M NaCl solution and 0.05 M NaCl solution as DS and FS, respectively, 0.34 kWh·m −3 of power can be produced in AB-PRO at R = 0.5, which is equivalent to that of a 120 m-high hydropower system.
Theoretical results suggest that AB-PRO is potentially superior to SS-PRO benefiting from the reversible mixing process and variable-pressure operation. The following sections will further assess its energy production when considering practical inefficiencies.

Impact of Inefficiencies
Inefficiencies as listed in Table 1, including pump efficiency (η P ), pressure exchanger efficiency (η PX ), energy recovery device efficiency (η ERD ), and pressure loss (P loss ) on the DS and FS sides, were considered when evaluating the practical performance of different PRO processes. The pressure drops on DS and FS sides were assumed to be the same (i.e., P loss = P F,loss = P D,loss ), and the P loss in AB-PRO was assumed to be 1/7 of that in SS-PRO as a shorter membrane module is allowed to be utilized in AB-PRO because of the smaller water recovery in each cycle. The parameters of η P of 0.85, η PX of 0.98, η ERD of 0.90, P loss of 0.2 bar in AB-PRO and P loss of 1.4 bar in SS-PRO were set as the baseline conditions. The impact of each parameter among η P , η ERD , and P loss was investigated by the control variate method.

Overall Impact of Inefficiencies
As revealed by Figure 3, inefficiencies induce significant reductions in the SEP and EPE of all three PRO processes. For SS-PRO with both flow modes, the energy production cannot compensate the unavoidable energy loss caused by device inefficiencies and pressure drops in practice at low R. An R of at least 0.1 is required to gain the net energy production as indicated by Figure 3a,b. When R ranges from 0.1 to 0.25, the SEP and EPE have an insignificant difference between the co-current SS-PRO and counter-current SS-PRO, which rises with the increase of R in both SS-PRO processes. At a higher R, the EPE of the cocurrent SS-PRO decreases due to the rapid decrease of ∆π at the outlet of pressure vessels while that of the counter-current SS-PRO remains stable (above 50%) benefitting from the more uniform distribution of ∆π along the membrane module.

Impact of Inefficiencies
Inefficiencies as listed in Table 1, including pump efficiency ( ), pressure exchanger efficiency ( ), energy recovery device efficiency ( ), and pressure loss ( ) on the DS and FS sides, were considered when evaluating the practical performance of different PRO processes. The pressure drops on DS and FS sides were assumed to be the same (i.e., = , = , ), and the in AB-PRO was assumed to be 1/7 of that in SS-PRO as a shorter membrane module is allowed to be utilized in AB-PRO because of the smaller water recovery in each cycle. The parameters of of 0.85, of 0.98, of 0.90, of 0.2 bar in AB-PRO and of 1.4 bar in SS-PRO were set as the baseline conditions. The impact of each parameter among , , and was investigated by the control variate method.

Overall Impact of Inefficiencies
As revealed by Figure 3, inefficiencies induce significant reductions in the SEP and EPE of all three PRO processes. For SS-PRO with both flow modes, the energy production cannot compensate the unavoidable energy loss caused by device inefficiencies and pressure drops in practice at low . An of at least 0.1 is required to gain the net energy production as indicated by Figure 3a,b. When ranges from 0.1 to 0.25, the SEP and EPE have an insignificant difference between the co-current SS-PRO and counter-current SS-PRO, which rises with the increase of in both SS-PRO processes. At a higher , the EPE of the co-current SS-PRO decreases due to the rapid decrease of ∆ at the outlet of pressure vessels while that of the counter-current SS-PRO remains stable (above 50%) benefitting from the more uniform distribution of ∆ along the membrane module. A significant energy loss of over 40% also occurs in AB-PRO when considering practical inefficiencies (baseline). As the is no longer approaching 0 in practical scenarios (e.g., the optimal is ~0.15 in the baseline case), the maximum hydraulic pressure that can be applied will be lower according to Equation (9), leading to more energy loss caused by under-pressurization (UP). Moreover, there is a larger salinity gradient between the recirculated solution and residual solution in both DS and FS tanks when increases, hence increasing the amount of entropy generation by mixing and energy loss. Furthermore, there is a tradeoff in AB-PRO: where a higher r leads to an increase of energy loss caused by UP and mixing, but a reduction of energy loss by pressure drops and device inefficiencies as the total volume of solution passing through the devices and membrane module is diminished (Equations (A8) and (A9) in Appendix B). Therefore, the SEP of AB-PRO was optimized by optimizing at each in this study, including the data in Figure  3c and other figures showing SEP of AB-PRO in the following sections. Although AB-PRO A significant energy loss of over 40% also occurs in AB-PRO when considering practical inefficiencies (baseline). As the r is no longer approaching 0 in practical scenarios (e.g., the optimal r is~0.15 in the baseline case), the maximum hydraulic pressure that can be applied will be lower according to Equation (9), leading to more energy loss caused by under-pressurization (UP). Moreover, there is a larger salinity gradient between the recirculated solution and residual solution in both DS and FS tanks when r increases, hence increasing the amount of entropy generation by mixing and energy loss. Furthermore, there is a tradeoff in AB-PRO: where a higher r leads to an increase of energy loss caused by UP and mixing, but a reduction of energy loss by pressure drops and device inefficiencies as the total volume of solution passing through the devices and membrane module is diminished (Equations (A8) and (A9) in Appendix B). Therefore, the SEP of AB-PRO was optimized by optimizing r at each R in this study, including the data in Figure 3c and other figures showing SEP of AB-PRO in the following sections. Although AB-PRO features a similar EPE (within a 5% difference) as the counter-current SS-PRO, its SEP and EPE are significantly higher than those of both SS-PRO processes at R < 0.4. In addition, the performance of AB-PRO is relatively stable. The change of EPE with R does not exceed 10%, which facilitates the flexible adjustment of R to meet the fluctuating electricity or pressure demand in application. When R = 0.5 in the baseline case, SEP and EPE are 0.19 kWh·m −3 and 56%, respectively. It should be noted that the value of SEP is related to the normalization method. The total energy production is normalized by the total volume of DS and FS in this study, but the SEP can also be defined as the energy production per volume of DS or FS to show the osmotic energy harvesting capacity if FS or DS is not limited. The results normalized by the initial DS volume ( Figure A1 in Appendix C) exhibit a higher SEP of 0.48 kWh·m −3 at R = 0.5 under the same conditions.

Impact of Pump Efficiency
By comparing the PRO performances at the same mode but different η P , the impact of pump on SEP and EPE is not significant as demonstrated in Figure 4. In all three PRO processes, an improvement of pump efficiency from 0.80 to 0.95 only results in an increase in EPE by less than 5%. The results are mainly contributed by the configuration of a highefficiency pressure exchanger (η PX = 0.98) which recovers most of the pressure generated by pumps and alleviates the negative effect of pump inefficiency, while the pump used in PRO is only for boosting the pressure after the PX, typically for several bars (<5 bar). The insignificant impact of η P implies more freedom in the choice of pump quality in practice. Moreover, as η P generally varies at various flow rates, the results also indicate a strong stability of the three PRO systems under different operating conditions. features a similar EPE (within a 5% difference) as the counter-current SS-PRO, its SEP and EPE are significantly higher than those of both SS-PRO processes at < 0.4. In addition, the performance of AB-PRO is relatively stable. The change of EPE with does not exceed 10%, which facilitates the flexible adjustment of to meet the fluctuating electricity or pressure demand in application. When = 0.5 in the baseline case, SEP and EPE are 0.19 kWh·m −3 and 56%, respectively. It should be noted that the value of SEP is related to the normalization method. The total energy production is normalized by the total volume of DS and FS in this study, but the SEP can also be defined as the energy production per volume of DS or FS to show the osmotic energy harvesting capacity if FS or DS is not limited. The results normalized by the initial DS volume ( Figure A1 in Appendix C) exhibit a higher SEP of 0.48 kWh·m −3 at = 0.5 under the same conditions.

Impact of Pump Efficiency
By comparing the PRO performances at the same mode but different , the impact of pump on SEP and EPE is not significant as demonstrated in Figure 4. In all three PRO processes, an improvement of pump efficiency from 0.80 to 0.95 only results in an increase in EPE by less than 5%. The results are mainly contributed by the configuration of a highefficiency pressure exchanger ( = 0.98) which recovers most of the pressure generated by pumps and alleviates the negative effect of pump inefficiency, while the pump used in PRO is only for boosting the pressure after the PX, typically for several bars (<5 bar). The insignificant impact of implies more freedom in the choice of pump quality in practice. Moreover, as generally varies at various flow rates, the results also indicate a strong stability of the three PRO systems under different operating conditions.

Impact of Energy Recovery Device Efficiency
As shown in Figure 1, an energy recovery device (ERD), such as a hydro-turbine (HT) or a pressure exchanger (PX), is installed in PRO processes to produce electricity or transfer the pressure to another system. The impact of on the SEP and EPE is significant, as demonstrated by Figure 5. With an increase in , the energy production performance is substantially improved in all three PRO processes, while the trends of performance curves for all the processes remain consistent with those in the baseline case. The maximum EPE increment in both SS-PRO and AB-PRO can approach 9% when increases from 0.85 to 0.95. By using a HT with a high efficiency of 0.95, totals of 0.19 kWh·m −3 , 0.22 kWh·m −3 , and 0.23 kWh·m −3 of electricity can be generated by the co-current SS-RO, counter-current SS-PRO, and AB-PRO, respectively, at DS = 1.2 M NaCl solution,

Impact of Energy Recovery Device Efficiency
As shown in Figure 1, an energy recovery device (ERD), such as a hydro-turbine (HT) or a pressure exchanger (PX), is installed in PRO processes to produce electricity or transfer the pressure to another system. The impact of η ERD on the SEP and EPE is significant, as demonstrated by Figure 5. With an increase in η ERD , the energy production performance is substantially improved in all three PRO processes, while the trends of performance curves for all the processes remain consistent with those in the baseline case. The maximum EPE increment in both SS-PRO and AB-PRO can approach 9% when η ERD increases from 0.85 to 0.95. By using a HT with a high efficiency of 0.95, totals of 0.19 kWh·m −3 , 0.22 kWh·m −3 , and 0.23 kWh·m −3 of electricity can be generated by the co-current SS-RO, counter-current SS-PRO, and AB-PRO, respectively, at DS = 1.2 M NaCl solution, FS = 0.05 M NaCl solution, and R = 0.60. Moreover, if there is another system requiring pressure, η ERD can be further enhanced to 0.98 by replacing the HT by a PX. In this case, the EPE of~59% can be maintained in both the AB-PRO and counter-current SS-PRO processes when R > 0.8. A higher EPE ranging from 60 to 66% can be realized in AB-PRO at R < 0.8, which is superior to that in the two SS-PRO processes, especially at low R. FS = 0.05 M NaCl solution, and = 0.60. Moreover, if there is another system requiring pressure, can be further enhanced to 0.98 by replacing the HT by a PX. In this case, the EPE of ~59% can be maintained in both the AB-PRO and counter-current SS-PRO processes when > 0.8. A higher EPE ranging from 60 to 66% can be realized in AB-PRO at < 0.8, which is superior to that in the two SS-PRO processes, especially at low .

Impact of Pressure Loss
Pressure loss ( ) caused by the friction of fluid in the flow channel (e.g., in membrane modules) of both DS and FS induces energy loss and lessens the net energy extractable by PRO. When assuming the same optimal hydrodynamic conditions (e.g., hydraulic diameter, flow velocity, and regimes) for SS-PRO and AB-PRO, the is proportional to the module length according to the Darcy-Weisbach equation [18,27]. In the following investigation, the of AB-PRO was assumed to be 1/7 of that of SS-PRO (Table 1) since batch PRO can use much shorter membrane modules to achieve a similar recovery in multiple-cycle operation. Moreover, the membrane module was assumed to have high mechanical strength, hence there was no membrane deformation in the PRO processes [38].
The results in Figure 6 indicate a significant decay of the energy harvesting performance when increases. In SS-PRO, a higher requires a higher threshold to gain net energy production and leads to lower SEP and EPE at the same . For example, if rises from 1.4 bar to 3.5 bar, the minimal for net energy generation increases from ~0.1 to ~0.2 in both SS-PRO processes. Meanwhile, the reduction of SEP ranges from 15% to 42% at different recoveries, and the highest EPE is only ~30% at = 0.5 in the cocurrent mode and ~40% at = 0.9 in the counter-current mode. On the contrary, the EPE can increase by 5-25% if can be reduced to 0.7 bar in SS-PRO through optimizing hydrodynamic conditions, feed spacer geometry, and so on [38,39]. The counter-current SS-PRO always outperforms the co-current SS-PRO at various and is capable of maintaining a relatively stable EPE of ~60% at > 0.30 and = 0.7 bar. At the of 0.5 bar, AB-PRO presents a better performance than SS-PRO (at the of 3.5 bar) and realizes an EPE of 40-50% benefitting from the shorter module length. However, alleviating from 0.2 bar (baseline) to 0.1 bar only results in a 2-5% increase in EPE in AB-PRO. Although the decrease of per-cycle water recovery ( ) from ~0.22 to ~0.15 because of the reduction in brings about less energy loss caused by UP and entropy generation, it also induces more energy loss owing to inefficiencies in AB-PRO, resulting in a slight improvement of EPE. When = 0.1 bar in AB-PRO and 0.7 bar in SS-PRO, the counter-current SS-PRO shows a small advantage compared to AB-PRO at > 0.3,

Impact of Pressure Loss
Pressure loss (P loss ) caused by the friction of fluid in the flow channel (e.g., in membrane modules) of both DS and FS induces energy loss and lessens the net energy extractable by PRO. When assuming the same optimal hydrodynamic conditions (e.g., hydraulic diameter, flow velocity, and regimes) for SS-PRO and AB-PRO, the P loss is proportional to the module length according to the Darcy-Weisbach equation [18,27]. In the following investigation, the P loss of AB-PRO was assumed to be 1/7 of that of SS-PRO (Table 1) since batch PRO can use much shorter membrane modules to achieve a similar recovery in multiple-cycle operation. Moreover, the membrane module was assumed to have high mechanical strength, hence there was no membrane deformation in the PRO processes [38].
The results in Figure 6 indicate a significant decay of the energy harvesting performance when P loss increases. In SS-PRO, a higher P loss requires a higher R threshold to gain net energy production and leads to lower SEP and EPE at the same R. For example, if P loss rises from 1.4 bar to 3.5 bar, the minimal R for net energy generation increases from~0.1 tõ 0.2 in both SS-PRO processes. Meanwhile, the reduction of SEP ranges from 15% to 42% at different recoveries, and the highest EPE is only~30% at R = 0.5 in the co-current mode and 40% at R = 0.9 in the counter-current mode. On the contrary, the EPE can increase by 5-25% if P loss can be reduced to 0.7 bar in SS-PRO through optimizing hydrodynamic conditions, feed spacer geometry, and so on [38,39]. The counter-current SS-PRO always outperforms the co-current SS-PRO at various P loss and is capable of maintaining a relatively stable EPE of~60% at R > 0.30 and P loss = 0.7 bar. At the P loss of 0.5 bar, AB-PRO presents a better performance than SS-PRO (at the P loss of 3.5 bar) and realizes an EPE of 40-50% benefitting from the shorter module length. However, alleviating P loss from 0.2 bar (baseline) to 0.1 bar only results in a 2-5% increase in EPE in AB-PRO. Although the decrease of per-cycle water recovery (r) from~0.22 to~0.15 because of the reduction in P loss brings about less energy loss caused by UP and entropy generation, it also induces more energy loss owing to inefficiencies in AB-PRO, resulting in a slight improvement of EPE. When P loss = 0.1 bar in AB-PRO and 0.7 bar in SS-PRO, the counter-current SS-PRO shows a small advantage compared to AB-PRO at R > 0.3, but the EPE difference between them does not exceed 3%. AB-PRO has an absolute predominance at R < 0.30. but the EPE difference between them does not exceed 3%. AB-PRO has an absolute predominance at < 0.30.

Conclusions
The performance of atmospheric batch PRO (AB-PRO) was explored and compared with conventional single-stage PRO (SS-PRO) with different flow schemes. Variable-pressure AB-PRO with an infinitesimal per-cycle water recovery ( ) can approach the thermodynamic maximum SEP and 100% of EPE under ideal conditions, while the efficiencies of two SS-PRO processes decrease with the increasing owing to the irreversible energy loss of under-pressurization caused by the constant-pressure operation. The impact of inefficiencies, including device deficiencies and pressure losses, was also investigated for all three PRO processes. In the practical case, a significant decay of performance was observed in all three PRO processes. Although counter-current SS-PRO shows a comparable performance with AB-PRO at high , the SEP and EPE of AB-PRO are significantly higher than that of SS-PRO at low . AB-PRO is capable of maintaining a relatively stable and high efficiency in the entire range of , which facilitates meeting the fluctuating energy or pressure demand in application by the flexible adjustment of . When utilizing a 1.2 M NaCl solution as the DS and a 0.05 M NaCl solution as the FS, the SEP and EPE of AB-PRO reach 0.19 kWh·m −3 and 56%, respectively, in the baseline case at = 0.5. For either AB-PRO or SS-PRO process, the advancement in pump and ERDs together with process design are critical to the further enhancement of energy production capacity. With a fixed high pressure exchanger efficiency of 0.98, pressure loss plays the most important role in the overall efficiency of PRO, followed by the ERD efficiency and pump efficiency.

Conclusions
The performance of atmospheric batch PRO (AB-PRO) was explored and compared with conventional single-stage PRO (SS-PRO) with different flow schemes. Variablepressure AB-PRO with an infinitesimal per-cycle water recovery (r) can approach the thermodynamic maximum SEP and 100% of EPE under ideal conditions, while the efficiencies of two SS-PRO processes decrease with the increasing R owing to the irreversible energy loss of under-pressurization caused by the constant-pressure operation. The impact of inefficiencies, including device deficiencies and pressure losses, was also investigated for all three PRO processes. In the practical case, a significant decay of performance was observed in all three PRO processes. Although counter-current SS-PRO shows a comparable performance with AB-PRO at high R, the SEP and EPE of AB-PRO are significantly higher than that of SS-PRO at low R. AB-PRO is capable of maintaining a relatively stable and high efficiency in the entire range of R, which facilitates meeting the fluctuating energy or pressure demand in application by the flexible adjustment of R. When utilizing a 1.2 M NaCl solution as the DS and a 0.05 M NaCl solution as the FS, the SEP and EPE of AB-PRO reach 0.19 kWh·m −3 and 56%, respectively, in the baseline case at R = 0.5. For either AB-PRO or SS-PRO process, the advancement in pump and ERDs together with process design are critical to the further enhancement of energy production capacity. With a fixed high pressure exchanger efficiency of 0.98, pressure loss plays the most important role in the overall efficiency of PRO, followed by the ERD efficiency and pump efficiency. SS-PRO is a continuous, constant-pressure process; hence the derivation of energy production is straightforward. The energy components in SS-PRO include:

Abbreviations
(1) Energy consumption by pumping (E pump ) (2) Energy recovered by the pressure exchanger (E PX ) (3) Energy consumption by compensating pressure loss (E loss ) (4) Energy production by an energy recovery device (E ERD ) The SEP is the net energy production normalized by the total volume of DS and FS.

Appendix B.2. AB-PRO
AB-PRO is a continuous but variable-pressure (P t ) process, hence the derivation of its energy production is based on the integration with respect to time (t). The duration time (t B ) of a batch of AB-PRO can be obtained by Equation (A7).
It was assumed that the total water permeation rate (Q P ) is constant during the AB-PRO process. Thus, regarding a constant recirculation flow rate of FS (Q F ), r is also constant.
The volume of solution passing through the devices and membrane module on the FS side (V pass,F ) and DS side (V pass,D ) in t B can be obtained by Equations (A8) and (A9), respectively.
The energy components in AB-PRO include: (1) Energy consumption by pumping (E pump ) (2) Energy recovered by pressure exchanger (E PX ) (3) Energy consumption by compensating pressure loss (E loss ) E loss = 1 η P P D,loss Q D t B + 1 η P P F,loss Q F t B (A12) (4) Energy production by energy recovery device (E ERD ) The SEP of AB-PRO can also be obtained by Equation (A6) combining with Equations (9)-(12) and further simplified to Equation (13) in the main text.
In the ideal case, not only η P = η PX = η ERD = 1 and P F,loss = P D,loss = 0 bar, but also r → 0 . The highest SEP of AB-PRO can be approached by the following equation, which is the same as the thermodynamic maximum SEP (Equation (2) The energy components in AB-PRO include: (1) Energy consumption by pumping ( ) (4) Energy production by energy recovery device ( ) The SEP of AB-PRO can also be obtained by Equation (A6) combining with Equations (9)- (12) and further simplified to Equation (13) in the main text.
In the ideal case, not only = = = 1 and , = , = 0 bar, but also → 0. The highest SEP of AB-PRO can be approached by the following equation, which is the same as the thermodynamic maximum SEP (Equation (2) in the main text).